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POLYNOMIAL CONSTRAINTS ON REPRESENTING 
ENTANGLED QUBITS AS MATRIX PRODUCT 
■ V STATES 

O . ANDREW CRITCH AND JASON MORTON 

>: 

q , Abstract. We quantify the representational power of matrix prod- 

uct states (MPS) for entangled qubit systems by giving polynomial 
expressions in a pure quantum state's amplitudes which hold if and 
only if the state is a translation invariant matrix product state or 
a limit of such states. For systems with few qubits, we give these 
equations explicitly, considering both periodic and open bound- 
ary conditions. Using the classical theory of trace varieties and 
trace algebras, we explain the relationship between MPS and hid- 
den Markov models and exploit this relationship to derive useful 
■ parameterizations of MPS. We make four conjectures on the iden- 

^ \ tifiability of MPS parameters. 

CN ; Matrix product states (MPS) provide a useful model of 1-D quan- 

^ ■ turn spin systems which approximate the ground states of gapped local 

Hamiltonians [16]. Accordingly the problem of classifying phases of 
matter for such chains has been reduced to understanding equivalence 
classes (such as under LU operations) in the space of quantum states 
representable as matrix product states [T71 [SJ 
CN ■ With periodic or open boundary conditions, we describe the closure 

of this space of states representable by translation invariant binary 
MPS as an algebraic variety. Our description is given as an ideal of 
polynomials in the state's amplitudes that vanish if and only if the state 
is a limit of MPS with N spins and D = d = 2 dimensional virtual and 
physical bonds. In small cases our description is complete. In Section 
[H we exhibit a polynomial which vanishes on a pure state if and only 
if it is a limit of binary translation invariant, periodic boundary MPS 
with N = 4, and a set of 30 polynomials which vanish when N = 5. 
We also obtain many linear equations which are satisfied for N up to 
12. In Section 121 Theorem [7] gives an analogous result for MPS with 
open boundary conditions and N = 3. Finally we examine cases where 
iV>0. 

Matrix product states bear a close relationship to probabilistic graph- 
ical models known as a hidden Markov models (HMM)[3]. In Section 
[3j we make this relationship precise by modifying the parametrization 
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of HMM to obtain MPS. We review the invariant theory of trace iden- 
tities and trace varieties that has been used to study HMM [I], and 
how these results apply to varieties of MPS. In particular we obtain a 
nice parametrization for translation invariant binary MPS with peri- 
odic boundary conditions. Finally in Section @] we suggest other such 
relationships between probabilistic graphical models and tensor net- 
work state models. Our results are complimentary to the connection 
between invariant theory and diagrammatic representations explored 
in [5] and the approaches to quantum state tomography for MPS de- 
veloped in [2113]. 

1. REPRESENTABILITY BY TRANSLATION INVARIANT MATRIX 

PRODUCT STATES 

First consider a translation-invariant matrix product state with peri- 
odic boundary conditions. Suppose the inner (virtual) bond dimension 
is D, the outer (physical) bond dimension is d, and there are N spins. 
Fix DxD complex parameter matrices Aq, . . . , Ad-i, defining the same 
DxDxd parameter tensor at each site. This defines the tensor network 
state, for ij 6 {0, . . . , d — 1}, 

(1) *= Yl te(A 1 ---A iN )\i 1 i 2 ---t N }. 

ii,...,i N 

Question 1. Fixing virtual and physical bond dimension, which states 
are matrix product states ? 

Including states which are limits of MPS, a precise answer to this 
question could be given as a constructive description of the set of poly- 
nomials / in the coefficients of \I/ such that /(^u,...^) = if and only 
if if> is a limit of MPS. This would describe the (closure of the) set of 
MPS as an algebraic variety. See [B] for background on varieties. 

Such a description is possible because of the way MPS are defined. 
Each coefficient ipi ly „^ N is a polynomial function of the parameters 
a rst in the D x D x d tensor A. Thus ([1]) defines a regular map 
# : C° 2d -> C dN , whose image we denote by PB(D, d, N), the set 
of tensors representable by translation-invariant matrix product states 
with periodic boundary conditions. Its closure PB(D,d, N) in either 
the Zariski or classical topology is an irreducible algebraic variety con- 
sisting of those tensors which can be approximated arbitrarily well by 
MPS. We can thus refine Question [1] as follows. 

Question 2. Fixing, D, d, and N, what polynomial relations must the 
coefficients of a matrix product state satisfy: what is the defining ideal 
ofPB(D,d,N)? 
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Figure 1. Translation-invariant MPS with periodic boundary. 

We primarily examine the fully binary case D = d = 2. The invari- 
ance of trace under cyclic permutations of the matrices A ix , . . . , A iN 
means we can immediately restrict to the subspace spanned by bi- 
nary necklaces (equivalence classes of binary strings under cyclic per- 
mutation). For iV = 3 physical legs, this is the coordinate subspace 
(V'ooo : ^100 : V'no : V'ni) an d all three-qubit states with cyclic sym- 
metry are matrix product states. For iV = 4 it is the six-dimensional 
coordinate subspace Ooooo : ^1000 : V'noo : V'ioio : V'mo : V'mi) and 
not all states are MPS (Theorem [3]). In the N = 5 case the 8 equiva- 
lence classes of coefficients under cyclic permutation are V'ooooo, ^ioooo> 
V'nooo, V'moo, ^niio, ^lini, V'ioioo, and Vnoio- For AT = 6, . . . , 15 the 




Figure 2. The eight binary necklaces for N = 5. 



dimensions of this "necklace space" are 14, 20, 36, 60, 108, 188, 352, 
632, 1182, and 2192. In general there are 



d-axy necklaces of length N, where (p is Euler's totient function. Thus 
translation invariant MPS with periodic boundary of length A" and 
physical bond dimension d live in a linear space isomorphic to C nd ^ N \ 
Naively we have 8 parameters in our 2x2x2 tensor A, but on each vir- 
tual bond we can apply a gauge transformation P{ ■ )P~ l for P G SL 2 
without changing the state[13j. Since SL 2 is 3-dimensional, we ex- 
pect PB(2,2,3) to be 5- dimensional. Counting this way, our expected 
dimension of PB(D,d, N) is mm{D 2 (d — 1) + l,nd(N)}. We expect 
PB(D, d, N) to be a hypersurface when this equals rid(N), which hap- 
pens first when (D, d, N) = (2, 2, 4). In this case our expectation holds: 

Theorem 3. A four-qubit state ^ is a limit of binary periodic transla- 
tion invariant MPS with N = 4 if and only if the following irreducible 
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polynomial vanishes: 

^1010^1100 - 2 ^noo - #1000^1010^100^1110 + ^looo^uooV'iiio 

- #iooo'0ioiq'0iiio + 2 ^oooo^?oio^iiio + IS^ooo^ioioV'iiooV'nio 

- #OOOoV'lOlQ / 011OoV , lllQ - 16 ^1000^1100^1110 + #0000^1010^1100^1110 

- #0000^1100^1110 - #0000^1000^1010^1110 + #0000^1000^1100^1110 

- ^oooo^nio + 2 ^i 2 ooo^?oio^iiii - Vtooo AicAiii - #1000^1010^1100^1111 
+ #iooqV>ioio'0iioo^iiii + 2^oooo^?oio^iioo^iin - #1000^1100^1111 

+ ^ooooV'iiooV'iiii - #1000^1010^1110^1111 + #0000^1000^010^1110^1111 
+ #iooqV'iiooV'iiioV ; iiii - #ooooV , ioooV'ioio?/'iiooV ; moV ; mi 

- 2V>ooooV'iooqV'iihAiii + 2^0000^1010^110^1111 ~ V'ioooV'iiii 
+ 2^0000 ^oooV'ioioV'mi - V'odooV'ioio^iiii- 

Hence, up to closure, the the set PB(2,2,4) of tensors that can be 
represented in the form ([T]) where A and A\ are arbitrary 2x2 ma- 
trices, is a sextic hypersurface in the space of 2x2x2x2 tensors 
invariant under cyclic permutations of the indices. The 30-term hyper- 
surface equation was found using a parametrization of the matrices that 
is similar to the birational parametrization of binary hidden Markov 
models given in pp. 

An example of a pure state on four qubits on which the polyno- 
mial / of Theorem [3] is nonvanishing, and so cannot be arbitrarily 
well approximated by such a matrix product state, is given by letting 

V'ioio = V'mo = - 1 / 4 and ^0000 = ^1000 = V'lioo = ^1111 = V 4 - In tnis 
example, /(^) = 2 -5 , which is the maximal value of /(^) attained on 
corners of the 6-D hypercube. 

The other cases with N < 15 when we expect PB to be a hypersur- 
face are when (D, d, N) = (2, 4, 6), (3, 3, 7), (5, 15, 12), (3, 71, 13), and 
(2, 296, 14). In general, we will need many more polynomials to define 
the space of matrix product states as their zero locus. As an example, 
consider PB(2, 2,5), which we expect to be a five- dimensional variety 
in the necklace space 

Theorem 4. Any homogeneous minimal generating set for the ideal of 
PB(2, 2, 5) must contain exactly 3 quartic and 27 sextic polynomials, 
possibly some higher degree polynomials, but none of degree 1, 2, 3, or 
5. 

Proof. Using the bi-grading of ITheorem 5\ we decompose the ideal / 
into vector spaces I r}S . For each (r, s) with ~(r+s) < 6, we select a large 
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number of parameter values A at random, and use Gaussian elimination 
to compute a basis for the vector space I r>s of polynomials vanishing at 
their images ^(A), which is certain to contain I rs . We then substitute 
indeterminate entries for A symbolically into the polynomials to ensure 
that they lie in I r>s , yielding a bihomogeneous basis for / in total degree 
< 6. □ 

This is interesting, because the variety only has codimension 3, but 
requires at least 30 equations to cut it out ideal-theoretically. Such a 
collection of 3 quartics and 27 quadrics was found and verified symbol- 
ically. Exact numerical tests (intersection with random hyperplanes) 
indicate that the top dimensional component of the ideal they gener- 
ate is reduced and irreducible of dimension 5, and is therefore equal to 
PB(2,2,5). 



1.1. Homogeneity and GL^-invariance. Note that the equation of 
ITheorem 31 is homogeneous of degree 6, and every monomial has the 
same total number of Is appearing in its subscripts. Every MPS variety 
will be homogeneous in such a grading: 

Proposition 5. For any D, d, N , the space of translation-invariant 
MPS limits with periodic boundary conditions is cut out by polynomials 
in which each monomial has the same total number ofOs, Is, . . . (d—l)s 
appearing in its subscripts. 

Proof. In fact we claim that the ideal of PB(D, d, N) is Z d -homogeneous 
with respect to d different Z-gradings degj for < i < d — 1, where 
degj(\l/j) is the number of occurrences of i in J. Since deg(ipj) : = 
- J^f^Q 1 degj(V'j) = 1, PB(D,d,N) is also homogeneous in the stan- 
dard grading. 

The usual parametrization where A , . . . , A^-i have generic en- 
tries, is Z d -homogeneous with respect to the grading above along with 
letting degj(A,-) = 1 when i = j and when i ^ j. Being a homoge- 
neous map, its kernel, the defining ideal of PB(D, d, N), is homogeneous 
in each of these gradings as well. □ 

In fact, the variety is homogeneous in a stronger sense because of an 
action of GLd on the parameter space of In the example above, the 
action is given by 

A (A , A x ) = (aA + bA u cAi + dAi) 
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which descends to an action on \l/ by 

V'ijfci = ^2 9i P 9jq9krgis^ 



9oo 9oi 
9w 9n 



pqrs 

pqrs 



The embedding (C*) d C GLd as diagonal matrices gives rise to the 
Z d homogeneity of the proposition above. 

1.2. Linear invariants and reflection symmetry. There are addi- 
tional symmetries peculiar to the case D = d = 2. For a generic pair 
of 2 x 2 matrices A ,Ai, there is a one-dimensional family of matri- 
ces P G SL 2 such that P~ 1 A i P are symmetric. Thus, a generic point 
^ e PB(2, 2, N) can be written as *(A ,Ai) with = Aj, and then 

*J = tr(n i6 J^i) = tr((n i6 jA-) T ) = tr (n j& eve r se(J)A') = ^reverse(J). 

This implies 

Proposition 6. // an N-qubit state ^ is a limit of binary periodic 
translation invariant matrix product states, then it has reflection sym- 
metry: ipj = VWrse(j) for all J. 

For N > 6, iV-bit strings can be equivalent under reflection but 
not cyclic permutation, so then PB(2,2,iV) admits additional linear 
invariants, i.e. linear polynomials vanishing on the model. For N = Q,7 
these are: 

PB(2,2,6): Vnoioo-V'i 10010 

PB(2, 2, 7) : Vmoioo - V'mooio and 0noiooo - ^noooio 

For small N we can find all the linear invariants of PB(2, 2, N) using 
the bigrading of lTheorem 51 as in the proof of ITheorem 41 Modulo the 
cyclic and reflection invariants, there are no further linear invariants 
for N < 7, but PB(2,2,8) has a single "non-trivial" linear invariant, 

"011010010 + "011001100 — "011001010 + 011101000 — "^11011000 — "011100100- 

For N = 9, 10, 11, and 12, PB(2, 2, N) admits 6, 17, 44, and 106 such 
non-trivial invariants, in each case unique up to change of basis on the 
vector space they generate. 



2. MPS WITH OPEN BOUNDARY CONDITIONS 

We now consider matrix product states with open boundary con- 
ditions, which are even more similar to hidden Markov models than 
the periodic version. Here the state is determined by two boundary 
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state vectors bo , b\ G C , along with the D x D parameter matrices 
A ,..., A d _ x of the MPS, by 

(2) *= J2 blA il ---A lN b 1 \i l i 2 ---i N ) 

ii,...,i N 

(3) = Yl te(BA il ---A iN )\i 1 i 2 ---i N ) 

ii,...,i N 

where B = b\b^ is a rank 1 matrix. We denote the set of states obtain- 
able in this way by OB(D, d, N), and its closure (Zariski or classical) by 
OB(D,d, N). We do not have the cyclic symmetries of the PB model 
here, so we consider OB(D,d, N) as a subvariety of C dN . If the Ai 
and fog have non- negative entries with row sums equal to 1, and b% is a 
vector of l's, then (J2J) is exactly the Baum formula for HMM, so in fact 
the model HMM(D, d, N) studied in [1] is contained in PB(£>, d, N). 

The expression ([3]) for \1/ is invariant under the action of SLp on 
the Ai and B by simultaneous conjugation. Thus, we may assume 
B is in Jordan normal form, i.e. a matrix of all zeroes except pos- 
sibly in the top left corner. As well, the map (B, A\, . . . , Ad) (->■ 
(t~ N B, tAi, . . . , tAd) preserves so discarding the case B = (which 
will not change OB) we can assume that the top left entry of B is 1. 
Thus \I> is determined by dD 2 parameters, the entries of the A4. In 
particular, OB(2, 2, 3) is parametrized by (a dominant map from) 8 
parameters, and lives in an 8-dimensional space. This parametrization 
still turns out still to be degenerate: 

Theorem 7. A three-qubit state ^ is a limit of N = 3 binary trans- 
lation invariant MPS with open boundary conditions if and only if the 
following 22-term quartic polynomial vanishes: 

V'oiiV'ioo ~ ^ooi^onV'ioo^ioi - V'oio^onV'ioo^ioi + ^OOO^OllV'lOl 
+ ^001^010^011^110 - ^000^011^110 - ^010^011^100^110 
+ ^ooiV'oio'0101'0110 + V'ooiV'iooV'ioiV'iio - V'oooV'ioi^iio - ^oWiio 
+ V'ooo^on^iio - V'ooiV'oioV'iii + ^oooV'oioV'oiiV'iii + ^oOlV'100'0111 
+ ^oioV'ioo'0111 - V'oooV'ohV'iooV'iii - V'ooiV'iooV'iii - ^oooV'ooi'0ioiV ; iii 
+ ^ooo^iooV'ioiV'iii + ^000^001^110^111 - ^000^010^110^111 

That is, the variety OB(2,2, 3) is a quartic hypersurface in C 8 cut out 
by the polynomial above. This polynomial previously appeared in the 
context of the HMM$F2\ . 

Proof. The map \I> and its image are homogeneous in the same grad- 
ing as described in ITheorem 5\ which we can use as in the proof of 
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ITheorem 41 to search for low degree polynomials vanishing on the vari- 
ety. When (D, d, N) = (2, 2, 3) the quartic from the theorem appears in 
this search. The quartic is prime, and therefore defines a 7-dimensional 
irreducible hypersurface in C 8 . On the other hand, the Jacobian of the 
map \I/ at a random point, e.g. the point where A , A\ have entries 
1, 2, 3, 4, 5, 6, 7, 8 in that order, has rank 7. Therefore OB(2, 2, 3) is of 
dimension at least 7, and contained in the quartic hypersurface above, 
so they must be equal. □ 

From [Theorem 7| we can derive conditions on OB(2, 2, N) for N > 4 
as well. There is a marginalization map from OB(2, 2, N) to OB(2, 2, 3) 
given by h-> J2\j\=n-3^ u f° r eacn ^ °f length 3, which commutes 
with the assignment b\ £) . g ■ Aj 3 ■ ■ ■ Aj N b%. In fact there are N — 
2 such marginalization maps, each given by choosing 3 consecutive 
indices I to marginalize to (summing over the remaining indices J). 
By composing these maps with the quartic polynomial above, we can 
obtain N — 2 quartic polynomials vanishing on OB(2, 2, N). 

By analogy to the case of hidden Markov models discussed in the 
next section, we make the following 

Conjecture 8. For N > 4, a generic N-qubit state can be recovered 
from its marginalization to any three consecutive states. That is, each 
marginalization map OB(2,2,iV) — > OB(2,2, 3) is a birational equiva- 
lence of varieties. 

The analogous statement with HMM in place of OB is shown to be 
true in [TJ. 

Conjecture 9. A generic N-qubit (D=d=2) translation invariant ma- 
trix product state ^ with open boundary conditions is determined up to 
phase by a reduced density operator which traces out all but a chain of 
three adjacent states, but no fewer. 

When the three adjacent states are qubits 1, 2, and 3 (the first three 
legs of the diagram), this amounts to saying that the group S l of unit- 
modulus complex numbers acts transitively on generic fibres of the 
real- algebraic map 



when restricted to OB(2,2,iV). Here the right hand side denotes an 
order 6 tensor with indices ji, j 2 , ja, k 2 , k 3 , and ^ denotes complex 
conjugation. 
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Figure 3. Parameterization of an MPS model as a com- 
plex HMM using complex E and T matrices with all row 
sums equal to z G C and copy dot (comultiplication) ten- 
sor (circle). Contraction of a region of the tensor network 
enclosed by a dashed line yields an A tensor. 



Conjecture 10. A generic N-qubit (D=d=2) periodic translation in- 
variant matrix product state is determined up to phase by a reduced 
density operator which traces out all but a chain of four adjacent states, 
but no fewer. 

Similarly, this amounts to saying that S 1 acts transitively on generic 
fibres of the map 



\i s ,...,i N 



fc2fc3/C4i5..ijv 



when restricted to PB(2, 2, N). 



3. Matrix product states as complex valued hidden 

Markov models 

We now explain how the polynomial in Theorem [3] was obtained, and 
connect the classical hidden Markov model and matrix product states 
through a reparametrizing rational map. The parametrization of the 
state \I/ is analogous to that of the moment tensor of a binary hidden 
Markov model used in [1] for symbolic computations. 

Let T be a 2 x 2 transition matrix and E a 2 x 2 emission matrix. For 
a (classical) hidden Markov model, T and E are nonnegative stochas- 
tic matrices (their rows sum to one), representing a four-dimensional 
parameter space. For PB, T and E will be complex with row sums all 
equal to some constant z G C, so they form parameter space isomorphic 
to C 5 . We parametrize the in terms of (T, E) by 



A = T, 



e i 






en 



^oo 

^10 



hi 



This is shown in Figure [3J grouping and contracting the E, T, and 
copy dot tensors into an A tensor yields a dense parameterization of 
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an MPS as depicted in Figure [U We then parameterize E and T with 
the five parameters u, v , b, c , z by setting 



Composing these formulae with the map (Aq,Ax) i-> ^ yields a re- 
stricted parametrization p^ : C 5 — > C 2 , whose image lies inside 
PB(2,2,iV). 

Proposition 11. The variety PB(2,2,A r ) is at most 5 -dimensional, 
and the image of our restricted parametrization pn is dense in it. 

Proof. Suppose \I> = ty(Ao,Ax) for Aq,A\ generic. First, we will trans- 
form the Ai by simultaneous conjugation with an element P of SL2 to 
a new pair of matrices A' , A[ such that A' Q has equal row sums and 
A[ = DA'q for a diagonal matrix D. Generically, Aq is invertible, and 
we can diagonalize the matrix AiAq , so we write U^AxA^U = D , 
and then U~ 1 AiU = DU^AqU. Next we find another diagonal ma- 
trix Di G SL 2 such that D^U -1 A U Di has equal row sums. Then 
let P = UDi and A\ = D^ 1 U~ 1 A i UDi, and we are done with our 
transformation. Now ^ = ty(A' Q ,A[) since simultaneous conjugation 
does not change trace products. But now letting z be the common 
row sum of A' , we can solve linearly for u, t>o, b, and Co to obtain 
* = p(u,v ,b,Co,z). □ 

In fact we know from exact computations in Macaulay2 |7J that the 
dimension dimPB(2, 2, N) = 5 for 4 < < 100. This is proven by 
checking that the Jacobian of p attains rank 5 at some point with 
randomly chosen integer coordinates, giving a lower bound of 5 on the 
dimension of its image. 

When parametrized using p, there are sufficiently few parameters 
and the entries of \1/ are sufficiently short expressions that Macaulay2 
is also able to compute the exact kernel of the parametrization, i.e. 
defining equations for the model. It is by this method that we obtain 
the hypersurface equation of lTheorem 31 as the only ideal generator for 



3.1. Identifying parameters of MPS. Determining the parameters 
of an MPS is related to quantum state tomography, and represents a 
quantum analog to the identifiability problem in statistics. The ex- 
tent to which the parameters can be identified can be addressed alge- 
braically. 



E 




T 




z + u — Co z — b + Co 
z — b — c z + b + c 




PB(2,2,4). 
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Given D x D matrices A± . . . Ad with indeterminate entries, we write 
Cd^ for the algebra of polynomial expressions in their entries that are 
invariant under simultaneous conjugation of the matrices by GL 2 . 

Sibirskii[15j, LeronpD], and Procesi[H] showed that the algebra Cr,,d 
is generated by the traces of products tr(A io ■ ■ ■ A in ) as n > varies. 
For this reason, Co,d is called a trace algebra. Its spectrum, Spec Cd^, 
is a trace variety. Since the coordinate ring of PB(.D, d, N) is a subring 
of C04, we have a map SpecCo^ — > C dN parametrizing a dense open 
subset of PB(D,d,N). 

In the case D = 2, Sibirskii showed further that the trace algebra C24 
is minimally generated by the elements tr(Ai) and tr(A^) for 1 < i < d, 
tr(AiAj) for 1 < i < j < d, and tr(AiAjA k ) for 1 < % < j < k < d. 

For d = 1 . . . 6, the number of such generators is 2, 5, 10, 18, 30, 47. 
In particular, when d = 2, the number of generators equals the tran- 
scendence degree of the ring, 5 = 8 — 3. This means SpecC 2)2 is 
isomorphic to C 5 , yielding for each N a dominant parametrization 
4>n : C 5 — > PB(2,2, N). Grobner bases for randomly chosen fibers in- 
dicate that for N = 4 ... 10, the map 4>n is generically fc-to-one, where 
k = 8, 5, 6, 7, 8, 9, 10, respectively. Continuing this sequence suggests 
the following. 

Conjecture 12. Using the trace parameterization <pN, for N > 5, al- 
most every periodic boundary MPS has exactly N choices of parameters 
that yield it. 

In other words, for N > 5, the parametrization (f)^ : C 5 ~ Spec 6*2,2 — > 
PB(2,2, N) is generically iV-to-1. Generically, the points of SpecC 2) 2 
are in bijection with the SX2-orbits of the tensors A. The conjec- 
ture implies that, up to the action of SL 2 , the parameters of a binary, 
D = d = 2 translation invariant matrix product state with periodic 
boundary are algebraically identifiable from its entries. 

4. Conclusion 

A conjectured dictionary between tensor network state models and 
classical probabilistic graphical models was presented in In this 
dictionary, matrix product states correspond to hidden Markov mod- 
els, the density matrix renormalization group (DMRG) algorithm to the 
forward-backward algorithm, tree tensor networks to general Markov 
models, projected entangled pair states (PEPS) to Markov or condi- 
tional random fields, and the multi-scale entanglement renormalization 
ansatz (MERA) loosely to deep belief networks. In this work we for- 
malize the first of these correspondences and use it to algebraically 
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characterize quantum states representable by MPS and study their 
identifiability. 

Acknowledgments. AC and JM were supported in part by DARPA 
under awards FA8650-10-C-7020 and N66001-10-1-4040 respectively 
We would like to thank J. Biamonte, J. Eisert, B. Sturmfels, F. Vac- 
carino, F. Verstraete, and G. Vidal for helpful discussions. 

References 

[1] Andrew J. Critch. Binary hidden Markov models and varieties. 

arXiv:1206.0500, 2012. 
[2] T. Baumgratz, D. Gross, M. Cramer, and M.B. Plenio. Scalable reconstruction 

of density matrices. arXiv preprint arXiv:1207.0358, 2012. 
[3] J. Biamonte, V. Bergholm, and M. Lanzagorta. Invariant theory for matrix 

product states. arXiv:1209.0631, 2012. 
[4] N. Bray and J. Morton. Equations defining hidden Markov models. In 

L. Pachter and B. Sturmfels, editors, Algebraic Statistics for Computational 

Bio., chapter 11, pages 237-249. Camb. Univ. Press, 2005. 
[5] X. Chen, Z.C. Gu, and X.G. Wen. Classification of gapped symmetric phases 

in one-dimensional spin systems. Physical Review B, 83(3):035107, 2011. 
[6] D Cox, J Little, and D O'Shea. Ideals, varieties, and algorithms. Springer- 

Verlag, New York, second edition, 1997. 
[7] Daniel R. Grayson and Michael E. Stillman. Macaulay2, a soft- 
ware system for research in algebraic geometry. Available at 

|http:// www.math.uiuc.edu/Macaulay2 / 
[8] M. B. Hastings. Entropy and entanglement in quantum ground states. Phys. 

Rev. B, 76:035114, 2007. 
[9] R. Hiibener, A. Mari, and J. Eisert. Wick's theorem for matrix product states. 

arXiv preprint arXiv: 1207. 6537, 2012. 
[10] U. Leron. Trace identities and polynomial identities of n x n matrices. J. Al- 
gebra, 42:369-377, 1976. 
[11] J. Morton. Tensor networks in algebraic geometry and statistics. Lecture at 

Networking Tensor Networks, Centro de Ciencias de Benasque Pedro Pascual, 

Bcnasque, Spain, May 2012. 
[12] L. Pachter and B. Sturmfels. Tropical geometry of statistical models. PNAS, 

101(46):16132-16137, 2004. 
[13] D. Perez-Garcia, F. Verstraete, MM Wolf, and JI Cirac. Matrix product state 

representations. Quantum Information & Computation, 7(5):401-430, 2007. 
[14] C. Procesi. The invariant theory of n X n matrices. Adv. Math., 19:306-381, 

1976. 

[15] K.S. Sibirskii. Algebraic invariants for a set of matrices. Siberian Mathematical 

Journal, 9:115-124, 1968. 
[16] F. Verstraete and J.I. Cirac. Matrix product states represent ground states 

faithfully. Physical Review B, 73(9):094423, 2006. 
[17] F. Verstraete, J.I. Cirac, JI Latorre, E. Rico, and M.M. Wolf. 

Renormalization-group transformations on quantum states. Physical review 

letters, 94(14):140601, 2005. 



